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Abstract 

We study the propagation of surface electromagnetic waves along the metallic surface covered by various layered 
dielectric structures. We show that strong radiative losses typical for the scattering of the surface wave can be consid- 
erably suppressed when single dielectric step is substituted by gradient index or periodic layered structure. 
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1. Introduction 

It is well known that the life time of surface elec- 
tromagnetic waves (surface plasmon polaritons - SPP) 
propagating along the metal-dielectric interface fTJ|2l[3] 
is rather short. Even in the case of ideal lossless metal, 
the propagation of the SPP through any inhomogeneity 
is accompanied with radiative losses 0. The simplest 
example is a transmission of the SPP through a single 
dielectric step, represented by the planar interface be- 
tween two dielectrics which cover the metallic surface 
J5]|6l. As was shown recently |7 8], SPP can radiate up 
to 40 per cent of its energy at the single interface. 

The origin of radiative losses lies in the different spa- 
tial distribution of the electric and magnetic field in the 
direction perpendicular to the metallic surface. To bal- 
ance the tangential components of E and H, the SPP 
radiates plane waves. 

The radiative losses can be suppressed for SPP prop- 
agating at the interface between metamaterials. The 
losses are almost negligible for the TE polarized SPP 
propagating along the surface of negative permeability 
material. The physical reason for this is that the spa- 
tial distribution of the TE polarized SPP depends only 
weakly on the permittivity €d of the dielectric J9). In 
Ref. Pol, anisotropic metamaterials were proposed to 
reduce radiative losses for the TM polarized SPP. 

The presence of plane waves makes the quantitative 
analysis of the propagation of SPP mathematically dif- 
ficult. Analytical calculation of the transmission and re- 
flection coefficients [6] of the SPP through discontinu- 



ities represents rather difficult task. Numerical methods 
0B are very useful for the quantitative analysis of the 
scattering. The problem is formulated for the scattering 
matrix of both SPP and radiative plane waves. With N w 
plane waves considered, the size of the scattering matrix 
increases to 2(1 + N w ) x 2(1 + N w ). 

In this paper, we analyze the propagation of the SPP 
along the metallic surface covered by various dielectric 
materials. Of special interest is the propagation through 
layered structure when the metal is covered by paral- 
lel dielectric strips of various with and permittivities. 
We consider gradient index structures, in which the per- 
mittivities in two adjacent strips differs in small value 
As, and periodic structures composed of alternating di- 
electrics. We calculate the transmission and reflection 
coefficients for the SPP and show that radiative losses 
can be considerably reduced when the single interface 
between two dielectrics is substituted by the layered 
structure in which the permittivity smoothly changes. 

The paper is organized as follows. In Section[2]we de- 
fine the model and introduce the parameters of the SPP. 
Section |3] discuss the the method of calculation of the 
transmission parameters for SPP propagating through 
the layered structure. In Section [4] we show that radia- 
tive losses can be dramatically reduced when the sin- 
gle interface ei/e2 is substituted by the gradient index 
structure (T3J, represented by N w planar layers with lin- 
early increasing dielectric permittivity. In Section [5] we 
discuss the propagation of the SPP through planar slab. 
Similarly to the case of the single interface, we find 
much higher transmission when the slab is modeled by 
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the gradient structure (sect. [6j. Propagation of the SPP 
through layered medium constructed by periodic alter- 
nation of two dielectrics e\ and ej. is analyzed in Section 
[7] We show that the propagation of SPP through such 
structure does not require an increase of radiative losses: 
radiation is the same, or even smaller, than that accom- 
panying the propagation through single interface. 

2. Definition of the modei 

Consider the metallic surface lying in the xy plane 
with metal filling the z < region. The metal is cov- 
ered by different dielectrics. All interfaces between two 
neighboring dielectrics are parallel to the yz plane. The 
SPP propagates along the metallic surface. The direc- 
tion of the propagation is perpendicular to dielectric in- 
terfaces. The electric and magnetic field of the SPP de- 
creases exponentially in the z direction as exp(-Kdz) for 
z > (dielectric) and exp(+/f m z) for z < (metal). The 
analytical form of parameters and K m is HI 



El. Here we only note that we will consider N w 
plane waves with the same frequency at but different z- 
component of the wave vector k z : 



k 

k ra = a, 

N w 



1,2,.. .N w . 



(5) 



Both N w and k zmdx represent the parameters of the 
model. Other components of the plane wave vectors are 
k y = and 



(6) 



Only waves with real k xa contribute to radiative losses. 
In what follows we consider the frequency of the SPP 



co = 0.23w p , 



(7) 



which corresponds to the frequency of visible light, and 
neglect absorption. Then, the metallic permittivity is 
e m = -17.9. Using Eq. [T] we express the wavelength A 
of the SPP: 



o 



Srf + E n 



and 



o 



(1) 



Here, k$ - a>/c and c is the light velocity, s,„ is a per- 
mittivity of metal, represented by lossless Drude for- 
mula, e m = 1 - (Dp/ co 2 and is dielectric constant of 
dielectric medium. The SPP propagates along the x di- 
rection. Then, its wave vector in the xy plane has only 
x-component, 

; 2 _ ; 2 s d s m 



Ed + £ m 

Equations ( T]2 1 give complete frequency dependence of 
the wave vector of the SPP. 

For the metal-dielectric interface only the TM polar- 
ized SPP can be excited [ 15 1. The intensity of the mag- 
netic field is 



h 



M)(0, 1,0) <D(x,y,f)e^ z 
M)(0, 1,0) ®(x,y,t)e +K '" z 



z>0, 
z<0, 



and the intensity of the electric field is 



Nq^ <D(x,y,f) (+iK d ,0,-k x )e *<> z /s d 
k 

M)r *(*.y.O (-iK m ,0,-k x )e +K " z /e m 
k 



(3) 

z>0, 
z < 0. 



Here, 0(x,y,f) = e*--^"") and z = yj/uo/eQ. The con- 
stant A^o normalizes the energy current in the x direction 

on. 

Similar formulas can be derived for the plane waves 
scattered at the metallic interface. Explicit expression 
for the electric and magnetic fields is given elsewhere 



lo p e d + e„ 
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where A p = 2nc/u) p » 94.2 nm is the wavelength cor- 
responding to the plasma frequency a) p /27V = 2 x 10 15 



3. Propagation of the SPP through planar interfaces 

Scattering of SPP on the single interface is accompa- 
nied by the radiation losses (8). Therefore, the scatter- 
ing matrix T is of the size of 2(N W + l)x2(N H + 1), where 
N is a number of plane waves, given by the z-component 
of the wave vector. The amplitudes of scattered waves 
are related by the scattering matrix 



S21 



S12 

S22 



(9) 



Here, A, and B, are the amplitudes of waves on the left 
and right side of the interface. Amplitudes A,- and B, 
belong to the plane waves propagating in the opposite 
direction. We have N w + 1 waves, (z = . . . N w ) with Ao 

(4) 

and Bq representing the amplitudes of the SPP. 

The scattering matrix expresses the waves propagat- 
ing away from the interface (B, A) through the ampli- 
tudes of the incident waves (A, B). Explicit form of 
scattering matrix is given in 0. The transmission 
and the reflection coefficients of the SPP incident from 
medium 1 are given by 



Ti^2 = |S 12 (00)| 2 , = |S 22 (00)| 2 , 



(10) 



and from medium 2 



so that 



T 2 ^ = |S 21 (00)| 2 , R 2 ^ 2 = |S„(00)| 2 



(11) 



Radiative losses into the first and second media are 

N, N 2 

S 1^1 = £ |S 22 (aO)| 2 , 5 ^ 2 - £ |S 12 C80)| 2 , (12) 



5 21 - '7~ 22 1 - r 22 /l(T 21 /lTi 2 T 22 /l + T 22 ) , (20) 

and the transmission coefficient of the SPP T = 
|S 2 i(0, 0)| 2 is expressed without any use of elements 
of A . Similarly, we express matrices 7~i 2 and S n = 

7125 21- 



S 2 ^ 2 = JjS„(aO)| 2 , S 2 ^ = J]|S 21 (/?0)| 2 , (13) S u = (t;^T 12 T 22 ^ + T' ]2 ) 



respectively. Here, Nj is a number of propagating plane 
waves (waves with real x component of the wave vector) 
in the ith medium. Conservation of energy requires 



l-»2 



1. 



(14) 



since the absorption is omitted. In this way we calcu- 
late transmission parameters of the surface wave prop- 
agating through the single planar interface. For more 
complicated structures which contain two or more in- 
terfaces, we apply the transfer matrix method. To every 
interface we assign it's transfer matrix defined as 



L BA 



S 12 ~ S 11^21 5 22 SllS^l 

-5 21 5 22 S~ x 



(15) 



The transfer matrix for the homogeneous dielectric 
spacing between two interfaces is 



A = 



A 
AT 



(16) 



Here A is a diagonal matrix containing phase factor 
£ ik m L^ w h ere £ j s me width of the dielectric slab. The 
transfer matrix of the entire structure containing n inter- 
faces is given as a product of successive transfer matri- 
ces 



T - T„„_iA„_iT, 



■ A2T21A1T10. 



(17) 



In the process of multiplication we have to keep in mind 
that some of the elements of the matrices Aj l correspond 
to waves that increases exponentially along the +x di- 
rection. Therefore the elements of T can become very 
large. To avoid numerical operations with A -1 we apply 
the normalization procedure lfT2l : consider the product 
T - T'AT. To calculate the transmission coefficient, 
we need the matrix 7~ 22 = which can be expressed 
as 



722 = (0 1) 



7-ii 
T21 



T22 



(18) 



This expression can be written in the form 



r 22 = (T 21 T 22 )(f 12T 2^)^T 22 (19) 



xOVTnT^ + Ty- 1 



(21) 



which determines the reflection coefficient R. This 
approach can be applied recursively to the multilayer 
structures. We calculate and Si" matrix ele- 

ments for (n + l)-th interface by 



o0>+l) 
13 21 
o(n+l) 
^ 11 



S^AY- 1 
AXY~ l A, 



(22) 



where 



X = T n Sf>+T n , Y = T 21 Sto+T22, (23) 

and Sjj , Sxv are scattering matrices for structure up to 
n-th interface. 



4. Gradient structures 

Previous scattering analysis [5, 6 7 8| showed that 
the propagation of the SPP through the single dielec- 
tric interface with high index contrast is accompanied 
by huge scattering losses. To suppress these losses we 
analyze layered structure which consists of N parallel 
layers of the same width. The permittivity of the j-th 
layer is 



£j — £a + As x i, 



(24) 



and the permittivity step As = (sb - £a)/N. This model 
mimics the gradient structure where the permittivity in- 
creases linearly from sa to sb- 

Figure [T] shows the transmission and the reflection 
coefficients and an amount of scattering losses for the 
normal incidence SPP propagating through such struc- 
ture for sa = 1 to Sb = 11.5. The structure is com- 
posed from N - 2Q dielectric layers of the same width. 
The transmission increases when the width of layers 
increases and saturates for sufficiently wide dielectric 
slabs. This effect is more significant for finer permit- 
tivity steps (Fig. [2|. It is obvious that radiative losses 
could be totally eliminated for perfectly smooth permit- 
tivity profile. We also see that the reflectivity of this 
structure is very small (Fig. [TJ. 




Figure 1 : The transmission coefficient T, the reflection coefficient R 
and radiative losses S for the normal incident the SPP propagating 
through the gradient structure with linearly increasing value of the 
permittivity. The permittivity changes from ea = 1 to eg = 11.5 with 
the increments Ae = 0.5 (spatial dependence of the permittivity is 
shown in the inset). The width L of each layer is measured relative to 
the wavelength A5 of the SPP. The transmission and scattering losses 
for the single interface £aI^b are T = 0.19 and 5 = 0.5, respectively. 

Of course, the limit of unity transmission coefficient 
can be reached only theoretically. With 20 permittiv- 
ity steps, considered in Fig. [T] the entire structure be- 
comes macroscopically large, when the width of layers 
is L — A5 ~ 245 nm. However, as shown in Fig. |2j con- 
siderable increase of the transmission coefficient can be 
obtained already with much thinner layers and/or with 
the use of only a few dielectric layers. 

5. Propagation of the SPP through dielectric slab 

Consider now the slab of dielectric material of width 
^siab with permittivity s B embedded into two semi- 
infinite layers with permittivity s A and study the prop- 
agation of the SPP across such structure. The transfer 
matrix (Tf) reduces to T — T AB A B T BA . 

Figure [3] shows the transmission of SPP through the 
slab with sa — 8 and eg = 1. The transmission and 
reflection coefficients exhibit typical Fabry-Perot res- 
onances which confirm that the transmission is deter- 
mined mostly by the ratio L s i a b/Ai of the slab width 



Figure 2: The transmission coefficient T for the normal incidence of 
the SPP propagating through the layered structure e, = ea + Ae X i for 
various permittivity increments Ae. Horizontal dashed line represents 
the transmission coefficient T = 0.19 for the single interface between 
media ea = 1 and eg = 11.5. 

to the wave length of the SPP. Similarly as in the case 
of the plane waves Ifl4l [TBI with maximal values for 
L = n x Ai/2 and maxima of the reflection R appear 
when L « (Ai/4) X (2n + 1). As shown in Fig. [3] the 
maxima of the transmission coefficients are accompa- 
nied with maxima of the scattering loses. 

The width dependence of the transmission parame- 
ters for the propagation of the SPP through the layer 
with higher permittivity, B — > A — > B, is shown in Fig. 
|4] Here, the transmission coefficient exhibits dramatic 
oscillations when the slab thickness varies. 

The difference between the two cases lies in the role 
of the plane waves. While in the first structure only a 
few plane waves are allowed to propagate through the 
slab (most of plane waves inside the slab are evanes- 
cent), the propagation of the SPP through the slab with 
higher permittivity is accompanied by a huge number of 
plane waves. The transmission coefficient of the plane 
wave a is determined by the product Lk xa . The radia- 
tive loses are enhanced each time when the x compo- 
nent of the wave vector of one of the plane waves fulfills 
the Fabry-Perot condition for the maximal transmission. 
With many plane waves with different value of k x inside 
the slab we observe huge oscillation of both the trans- 




Figure 3: The transmission parameters of the normal incidence of 
the SPP propagating through the dielectric slab with the permittivity 
eg = 1 sandwiched between two dielectric media with the permittivity 
ea = 8. The thickness of the slab is shown relative to wavelength of 
the SPP Aj . Fabry-Perot resonances of both the transmission and the 
reflection are clearly visible. Note that maxima of the transmission 
coefficient are accompanied with maxima of the radiative losses. 



mission coefficient of the SPP and of radiative losses S 
(not shown in Fig. [4]). Smooth L s i a b - dependence of the 
transmission coefficient is observed only for very thin 
slabs 



6. Quadratic structures 

The propagation of the SPP through the slab, dis- 
cussed in the previous section, can be enhanced with 
the use of gradient structures. As an example, we con- 
sider layered structure consisting of N dielectric layers 
with permittivities 

St = 1 + (2.75 - //4) 2 , i = 1 ... AT = 22. (25) 

The permittivity varies from the maximal value 8.56 to 
minimal value 1 in the middle of the structure (inset of 
Fig. [5j. The width L of each slab is again considered 
relative to the length of the SPP for the dielectric per- 
mittivity s d = 1 (A; = 632 nm). 



Figure 4: The transmission coefficient T, the reflection coefficient R 
and scattering losses for the propagation of the SPP through the dielec- 
tric slab with the permittivity eg = 8 sandwiched between dielectric 
media with the permittivities ea = 1 ■ The thickness of the is shown 
relative to wavelength of the SPP Ag . 



The transmission parameters of the SPP are shown in 
Fig. 5 Contrary to the propagation through the slab 
(Fig. 3 i, no Fabry-Perot oscillations are observable. In- 
stead, we see that the transmission coefficient increases 
continuously and the reflection coefficient is almost neg- 
ligible. 

The inverse structure, defined by the quadratic in- 
crease of permittivity 



Si = 8.5625 - (2.75 - i/Sy i = ... 44, 



(26) 



1 Note that the wave length As of the SPP is considerably smaller 
than litjkx for any accessible plane wave. 



exhibits oscillations of the transmission coefficient, 
which are, however, much more moderate than that 
shown in Fig. [4] for a single slab. The number of os- 
cillations is proportional to a number of permittivity 
steps. The transmission coefficient and radiation losses 
are anti -correlated: small radiation losses correspond to 
large transmission coefficient and vice versa. 

7. Propagation of the SPP through periodic struc- 
tures 

In this section we consider the propagation of the SPP 
through the layered structure built by N periods of lay- 
ers ABAB . . .AB with permittivities s A and s B . The pe- 



T | i | I | I | i | I | r 




L/A 



Figure 5: The transmission coefficient 7\the reflection coefficient 
R and radiative losses S for the normal incidence SPP propagating 
through the layered structure with quadratic change of the dielectric 
permittivity. Permittivity ej = 1 +(2.75-i'/4) 2 for i = ... 22 is shown 
in the Inset. The width of all layers is same and is shown relative to 
Al. 

riodic structure can be either sandwiched between two 
semi-infinite media A or between medium A and B form 
the left (right) hand side. 

We find that the SPP propagating through such struc- 
ture can radiate less energy than the process of the prop- 
agation through a single interface AB. Also, radiation 
losses do not increase as the number of stack increases. 
Waves radiated at the first interface are used in the pro- 
cess of scattering at next interface. 

Typical transmission properties of the periodic stack 
are shown in Figs. [7] and [8] The structure consists of 
N - 30 periods AB. The width of each layer is propor- 
tional to the wavelength of the SPP: 



and permittivities ea = 1 and sb - 2. 

Figure [7] presents the transmission parameters for 
the propagation of the SPP through N - 30 periods 
AB. Both left and right medium have the permittivity 
sa = 1. Transmission and reflection coefficients exhibit 
typical oscillations when the layer widths increase. For 
La/Aa > 0.22, we found the gap in which no propaga- 
tion of the SPP is possible. Similarly to plane waves, the 




L/A, 



Figure 6: The transmission coefficient T, the reflection coefficient 
R and scattering losses S for the normal incidence SPP propagating 
through the layered structure with quadratic change of the dielectric 
permittivity s d = 8.5625 - (2.75 - i/8) 2 for i = . . . 44. Width of all 
layers is same and is shown relative to As . 

origin of the transmission gap lies in the wave character 
of the SPP. However, there are two differences between 
the gap for SPP and for plane waves in one dimensional 
photonic crystals. First, the reflection coefficient of the 
SPP does not reach the unity inside the gap. A small 
radiative losses are still observed since there is no total 
gap for the TM polarized plane waves IT6l . Second, 
since plane waves must assist the propagation of the 
SPP through each interface inside the layered structure, 
we expect that the transmission coefficient of the SPP is 
considerably reduced for layers widths where the prop- 
agation of majority of plane waves is restricted. Also, 
the position of the gap for the SPP does not scale with 
the frequency because of the non-linear frequency de- 
pendence of the wave length A. 

Figure [8] shows similar data for the propagation 
through N — 30 periods AB but with right medium 
sb = 2. The radiation losses are small for any width 
of layers. Plane waves radiated at the first interface into 
the material B assist the transmission of the SPP through 
next interfaces BA so no additional waves should be ra- 
diated. Also, for thin layers, the transmission coefficient 
is higher than the transmission coefficient for the single 
interface AB (T = 0.85). 




Figure 7: The transmission coefficient T, the reflection coefficient R 
and scattering losses S for the normal incidence of the SPP propagat- 
ing through the periodic structure of two alternating media as a func- 
tion of layer thickness. Layers are arranged in order ABAB . . .ABA, 
number of periods is N = 30. Permittivities are £,4 = 1 and Eg = 2, 
N = 30. Widths of layers are given by Eq. \27) . 



The last two Figures , [9] and [T0| demonstrate the trans- 
mission parameters for the SPP propagating through pe- 
riodic layered structure with single impurity represented 
by additional layer B or A in the middle of the structure. 
Similarly to the plane wave photonic crystals, impurity 
creates isolated level in the band gap, where the reflec- 
tion coefficient decreases. Figure [9] shows that this de- 
crease of the reflection should not automatically imply 
a good transmission of the SPP; instead, decrease of the 
reflection is accompanied with plane wave radiation and 
the transmission coefficient is small. The transmission 
of SPP is better for the A impurity layer included into 
the layered media with slower permittivity contrast. For 
instance, Fig. 10 shows well pronounced impurity in- 
duced transmission peak (T as 0.7 inside the gap for 
sa/sb = 3/5 with impurity A located in the middle of 



the periodic structure. 



8. Conclusion 



We analyzed the propagation of the surface plasmon 
polariton - SPP - through various layered structures in 
which the metallic surface is covered by strips of dielec- 



Figure 8: The transmission coefficient T, the reflection coefficient R 
and scattering losses S for the normal incidence of the SPP propagat- 
ing through the periodic structure of two alternating media as a func- 
tion of layer thickness. Layers are arranged in order ABAB . . .AB, 
number of periods is N = 30. Permittivities are ea = 1 and Eg = 2. 
Widths of layers are given by Eq. \21\ . Note that the transmission 
through such structure can be larger than the transmission through 
single interface, Tq = 0.86. 



trie materials with different dielectric permittivity. We 
showed that the transmission coefficient of SPP can be 
considerably increased when the single permittivity step 
with high permittivity contrast is substituted by a gradi- 
ent structure, represented by N thin dielectric slabs with 
increasing dielectric permittivity. 

We showed that the propagation of the SPP through 
more than one dielectric interface must not necessary re- 
quire an increase of the radiative losses. Plane waves ra- 
diated in one interface assist in the process of the trans- 
mission of the SPP through the next interface. 

Although no absorption is considered in our analysis, 
we expect that qualitatively the same results in absorb- 
ing structures. Absorption only reduces the transmis- 
sion coefficient for the SPP and radiative losses. 

This work was supported by Project VEGA 0633/09. 
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